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Abstract. We discuss the property of the number density of a fluid of particles 
living in a curved surface without boundaries to be constant in the thermodynamic 
limit. In particular we find a sufficient condition for the density to be constant along 
the Killing vector field generating a given isometry of the surface and the relevant 
necessary condition. We reinterpret the effect of a curvature on the fluid in a physical 
way as responsible of an external "force" acting on the particles. 
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1. Introduction 

The physics of fluids of particles living in surfaces is a well known chapter of surface 
physics. A special role is played by low dimensional exactly analytically solvable fluids, 
as they inform approximate solutions in higher dimensions and general sum rules. In 
the statistical mechanics of continuous fluids, those where the particles are allowed to 
move in a continuous space, one finds exact solutions for various Coulomb fluids pp. 
For example the one-component Coulomb plasma (OCP) is exactly solvable in one- 
dimension [2]. In two dimensions B. Jancovici and A. Alastuey [31 H] proved that the 
OCP is exactly solvable analytically at a special value of the coupling constant, in their 
1981 work. Since then, a growing interest in two-dimensional plasmas has lead to study 
this system on various flat geometries El [7] and two-dimensional curved surfaces like 
the cylinder [9], the sphere [El EU H21 EH QI] , the pseudosphere pU M\ E] , and the 
Flamm's paraboloid [T8~l [19]. Among these surfaces only the last one is of non-constant 
curvature. The statistical mechanics of liquids and fluids in curved spaces is a field of 
growing interest [20] . 

Here we do not restrict ourselves to those exactly solvable cases but want to find 
a general property of any given fluid living on a curved surface without boundaries. A 
homogeneous fluid living in a plane (or in general an Euclidean space) is known [21] to 
have a constant density. This same conclusion holds for a (non-ideal) fluid living in a 
surface of constant curvature in its thermodynamic limit 0- In this paper we will state 
what can be said about the constancy of the density for a fluid living in a Riemannian 
surface without boundaries and embeddable in the three dimensional Euclidean space, 
in its thermodynamic limit. It is obvious that an ideal fluid (a gas) has a constant 
density on any surface and whether or not we are in the thermodynamic limit. But 
what can be said about the non-ideal fluid? 

The study of Ref. [18] showed that the OCP in the Flamm's paraboloid is indeed 
homogeneous. We expect this occurrence to be due to the long range nature of the 
Coulomb potential and argue that it cannot hold in general for other choices of the pair 
potential or of the surface. 

In this work we will give a physical interpretation to the curvature of the surface 
as an external "force" guiding the particles of the corresponding "flat" fluid. We will 
show that the Coulomb potential has to be a function of the geodesic distance between 
the charges and we will restrict to a definition of a fluid as one made of particles with 
a pair interaction potential which is a function of the geodesic distance between the 
two particles. We will then find a necessary and sufficient condition for the density 
multiplied by the square root of the determinant of the metric tensor to be constant 
along a certain direction. We will show how this condition holds true both for non- 
quantum and quantum fluids. 

The paper is organized as follows: in Sec. [2] we state the problem we want to solve 
at the level of the non-quantum fluids; in Sec. |3]we reformulate the problem in such way 

| The notion of thermodynamic limit will become clear further on in the paper. 
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as to make explicit the physical interpretation of the curvature of the surface; Sec. H]is 
devoted to the quantum fluid formulation of the problem; Sec. [5] is for final remarks. 

2. Statement of the problem 

Given a non-quantum fluid of point wise particles living in a surface S embeddable in 
the three dimensional Euclidean space (note that we will not take under consideration 
those surfaces deriving from a Riemannian metric but not embeddable and those not 
deriving from a metric) and without boundaries one can define the canonical ensemble 
particle number density as [2T] 



where N is the number of particles confined in the region Q, (3 = 1/ksT with 
Boltzmann's constant and T the absolute temperature. The potential energy of the fluid 
is V = J2i<i<j<N v {d( ( iii Qj)) where v is the pair-potential and d(q, q') is the geodesic 
distance between the two points q and q'. The surface is defined by a metric tensor g a p 
so that the square of the proper length of the infinitesimal line element is given, using the 
usual Einstein's summation convention, by ds 2 = g a/ 3(q)dq a ®dqP where ® is the usual 
tensor product. We denote with g(q) = det 1 \g a p((i) 1 1 the Jacobian of the transformation 
from a locally flat reference frame to the local coordinates system on the surface. Here 
we use a coordinate basis {e a = d qa } so that q = q a e a and the symbol d stands for the 
exterior derivative. As usual we use upstairs Greek indexes for contravariant components 
and downstairs Greek indexes for covariant components, and we use a downstairs roman 
index to denote the (distinguishable) particle number. The symbol A indicates the usual 
wedge product. In the following we will call vol(f2) = J n \J g(q) A„ =1 dq a the volume of 
the region Q. 

The problem we want to discuss is the one of finding continuous transformations 
that leave unchanged the density p(q) in the thermodynamic limit. Here we think of 
the surface S as an embeddable one without boundaries. And by thermodynamic limit 
we mean that if S extends to infinity, vol(f2) — > oo with p = iV/vol(f2) kept constant or 
if S is closed, Q — > S with p = N/vol(S). We want to answer the question: "when is 
p(q) constant on S in the thermodynamic limit?". 

The number density satisfies the following normalization condition 




(1) 



Z 




i=l 





(3) 



So when the density is constant in the surface we must have p = p. 
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3. Reinterpretation of the curvature 

Choosing the coordinate basis so that £ = d q a is a Killing vector field [22] generating an 
isometry, then g, a = 0, where we use the usual comma convention to indicate a partial 
directional derivative. We know that if p is the momentum of a free particle on S then 
p ■ £ is a constant of motion p a (p ■£)' a = 0, where we use the usual semicolon convention 
to indicate a covariant derivative. The ideal gas has constant density on every surface 
regardless of the curvature and of the thermodynamic limit. We thus have to worry 
about the term exp(— f3V). Now, if one moves the N particles at qi, . . . , along the 
vector field £ the geodesic distances among the system of particles will stay constant 
as well as the potential energy V. We then have proven that given a Killing vector 
field d q a then p a = 0. Strictly speaking before taking the thermodynamic limit, the 
domain has boundaries, and close to these, one might not be able to move the particles 
along the Killing vector field, invalidating the conclusion near the boundary. When 
taking the thermodynamic limit, one needs to be able to quantify if these boundary 
effects will be negligible or not, and how deep they can affect the bulk of the system. 
This depends of the pair potential v and on the surface. In the flat space it is well 
known that the boundary effects are negligible (for suitable short-ranged potentials and 
for the Coulomb potential for globally neutral systems to have screening). But for a 
general curved surface, a proper study of what happens in thermodynamic limit with 
this boundary effect is needed and it will certainly impose additional conditions on the 
pair potential v, and probably also on the surface, to keep valid the conclusion that 
p >ol — 0. The conditions on the surface might appear for example in cases similar to the 
pseudosphere, where it has been shown that boundary effects can be of the same order 
of magnitude as the bulk properties (see Refs. fT^ fT6| fT7]). So, additional work in this 
direction is needed. 

This is clearly only a sufficient condition but it is enough to say that on the sphere 
(or the plane), the surface of constant curvature [23], where £ = d^, with ip the azimuthal 
angle, the density will be constant in the thermodynamic limit. One, in fact, has that 
the density is constant along parallels. And this, given the symmetries of the sphere, 
means that the density is indeed everywhere constant over the whole sphere, with p = p. 

On the other hand a necessary condition can be expressed as follows: Say that we 
find a coordinate system such that, for all v , {yfgp), a = then in particular for v = we 
have p =constant and g >a = 0. For the Flamm's paraboloid [18] we can say that there 
certainly exists a fluid (at least one v) such that (yfgp), r 7^ since d r is not a Killing 
vector of the surface and g iT ^ 0. And we know [18] that the OCP is an example. 

The problem then reduces to understand what can be said about surfaces of non- 
constant curvature. Note that we can as well rewrite Eq. ([Q) as follows 

, r P ^[v(qi,,q i v)+E, , Ii«q.;ffl]n Af A 2 da a * 

v^rtu) = N T em, q Wi¥ J' A ' (4) 

where 0(q; 0) = — [lng(q)]/2/3 is an "external potential". A form which suggests, on 
physical grounds, a local dependence of the density on the curvature. The fluid is seen 
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in this formulation as living on a "flat space" , the two dimensional space determined by 
the local coordinates chart {q l ,q 2 ) used in the surface, subject to an external potential 
induced by the metric. This suggestive reinterpretation of the problem can sometimes 
lead to a wrong intuition. For example we know that the OCP on the Flamm's 
paraboloid (see Sec. 4.2.4 of Ref. [IB]) has a density that is everywhere constant 
even if this surface is only asymptotically flat but curved near the "horizon" , the scalar 
curvature being proportional to the Euclidean distance r from the origin to the power 
of minus three. Whereas the constancy of the density along the azimuthal direction 
(p has to be expected from the sufficient condition stated above, the constancy of the 
density along the radial r direction is not at all intuitive, even more so at the light of 
the discussion which follows. 

For a surface with a conformal metric g a p = \J g(<l)o~ai3, l§l the scalar curvature R 
can be written as 

R( q ) = e^W/?A flat 0(q) , (5) 

where Afl at = + d 2 2 is the flat Laplace's operator. The external "force" acting on the 
particles due to the curvature is then —R exp(— 0(f)) / (3. For the Flamm's paraboloid [18] 
the force acting on the charges turns out to be 4/[/3s(l + s) 2 ] where s = \J (q 1 ) 2 + (q 2 ) 2 - 
As we already mentioned above, in this case, the OCP shows a constant density in the 
surface. In Section 13.21 we show that in general it has to be expected a non-constant 
density. 

On the other hand the formulation of Eq. (j3J) suggests that certainly ^fgp is a more 
fundamental quantity than just p itself to look upon. 



3.1. The Coulomb pair potential 

Here we want to show that the Coulomb potential between two charged particles living 
in a given surface S has to be a function of the geodesic distance between the charges 

m n m mm m- 

The Coulomb potential is defined by the Poisson's equation, 

A q ^(q,q') = -2vr5 (2) (q,q') , (6) 

where A q is the Laplace-Beltrami operator and <5^(q, q') = 5^(d(q, q')) the Dirac delta 
function, in the surface S. The Laplace-Beltrami operator is invariant to isometries. 
This means that if the charge at q and the one at q' are moved along the vector field of 
an isometry the Laplace-Beltrami operator will not change. Neglecting eventual additive 
functions which have a null Laplacian we must have 

vcoul = /(d(q, q')) ■ (7) 

§ Note that the following are all surfaces of this kind: the sphere embedded in three dimensional 
Euclidean space ^fg = 4/(1 + s 2 ) 2 , the pseudosphere embedded in three dimensional Minkowski space 
y/g = 4/(1 — s 2 ) 2 , the cylinder embedded in three dimensional Euclidean space y/g = 1, and the Flamm's 
paraboloid embedded in three dimensional Euclidean space ^fg = (1 + 1/s) 4 . Here s = \J {q 1 ) 2 + (q 2 ) 2 - 
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For example on the sphere pTO] of radius R one finds f(x) = — \n(2Rsm(x/2R)/L) with 
L a length scale. The conclusion of Eq. (J7J) is in agreement with Fermat's principle for 
light propagation [21]. 

3.2. The Coulomb fluid 

For an open surface with a conformal metric g a p = (^/g(s)/s)5 a p, s 6 [0, +oo[ the 
Laplace-Beltrami operator can be rewritten as 

Af = -^W , (8) 

where A flat is the usual Laplace operator in flat space (x = s cos <p,y = s sirup). We can 
then introduce a complex coordinate z = se %Lp and the Laplacian Green's function ([6]) 

A flat fcw(0, <p), Oo<A>)) = -2tt-S(s - s )S(<p - p ) (9) 

can be solved as usual, by using the decomposition as a Fourier series. Since (J5]) reduces 
to the flat Laplacian Green's function, the solution is the standard one 

vcoui((s,<p), Oo<A))) = ^2 ~ ( ~) cos K<P - ¥>o)] + v (s, s ) , (10) 

n=l n ^ S> ^ 

where s> = max(s,s ) and s< = min(s, s ). The Fourier coefficient for n = has the 
form 

/ x J a,]" Ins + b£ s > s 
Vo(s,s ) = < _ _ , 11 

a In s + b s < s 

and it has to satisfy the boundary conditions that vq should be continuous at s = so, 
Oq In s + &o = % In s + , and its derivative discontinuous due to the Dirac's delta in 
(EJ), (Iq / so — Oq / so = — l/s . Summing explicitly the Fourier series ( ITUjl and requiring 
additionally that the Coulomb potential Vcmu{ s u s z) ^ e symmetric under exchange of 1 
and 2 we find 

v Co ui(s,ip; s ,fo) = -In 77 — + « , (12) 

with /i(s, So) = 1 or h(s, Sq) = y/sso and a a constant. Here if we imagine the plasma 
confined into a disk of radius R we can choose 

vcoui{s,ip; s , <A)) = -hi j-, — ^\ + b, (13) 

with h(s, s ) = -R and b = a — In R or /i(s, so) = ^sso and 6 = a so that if we rescale 
all the s into As and R into Ai? the Coulomb potential does not change apart from an 
additive constant. Imagine now we are on a plane [3], then h(s, sq) = R. Then in the 
definition of the density (JT]) at any temperature we can change integration variables in 
the numerator from (sj, (fi) to (xi = Sie % ^ Vt ~ Vl \ yi = (pi — cpi) for i = 2, 3, . . . , N with 
Jacobian 1. Calling Vb = Vb(s/R) = p J n ^ Vcoui(s, ip; s', tp') \J g(s') ds'ip' the neutralizing 
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background potential and vq the self energy of the background we can write 



R 



x 



R i>j>2 k=2 
e -Pv b { Xk e-**/R) ^ g ( Xke -iy k)dxkdyk . (14) 

The integral does not depend on tpi so p(s\, tfi) = p{si). Now we can make a change of 
variables where s^ — > s^/si for k = 2, 3, . . . , N and Rj s\ — > T so that 

N 

I \Xk - 1| \' " 

X 

T i>j>2 k=2 



e"^ (Sfe/T) v^I^O "^fcdj/fc . (15) 

On a plane ^/fl'(ssi) = ssi so that in Eq. (fT5l) there is a multiplicative factor s^ -1 . So 
in the thermodynamic limit T — oo and JV-j>oowe can say that p(si) =constant since 
we know that we must have a well defined thermodynamic limit. The same conclusion 
holds on a pseudosphere (see Sec. 4.3.2 of Ref. [IB]), on a cylinder (see Eq. (12a) of 
Ref. pj), and on a Flamm's paraboloid (see Sec. 4.2.4 of Ref. [IB])- In these cases the 
explicit analytic expression of the density has been determined for the finite system as 
a function of the properties of the surface at the special value of the coupling constant 
fiq 2 = 2. To the best of our knowledge there aren' t any analytical results, in the 
literature, where the OCP has been found to have a non-constant number density in the 
thermodynamic limit on a given curved surface, probably one has to resort to numerical 
simulations [25]. It certainly has to be expected that in a general curved surface the 
OCP in the thermodynamic limit may have a non-constant density otherwise it would 
mean that an OCP in the plane has a uniform density for an arbitrary external field. 
It might actually be true that the effects of the metric and the background potential 
cancel one another when the potential is determined by Poisson's equation, but if it's 
true, it will be necessary to solve for the potential in more detail to prove it. 

4. The quantum case 

For the quantum fluid we find for the canonical ensemble distinguishable density matrix 
(the full density matrix for a system of Bosons or Fermions is then obtained by 
symmetrization or anti-symmetrization respectively) [26] 

Pd(Q', Q; P) = J Pd(Q', Q((M - l)r); r) • • • p D (Q(r), Q; r) x 

M-l 2N 
j=l a =l 

where as usual we discretize the imaginary time in bits r = hfi/M and Q = (qi, . . . , q^v) 
with 

9(i) = det HiMQ( ir ))ll ' ( 17 ) 
<MQ) = &*ift(qi) ® • • • ® 9a N /3 N (qN) , (18) 
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to get to the path integral formulation and in the small r limit for particles of unitary 
mass follows 

p(Q(2r),Q(r);r) = (2irh)- N g^y/D(Q(2T), Q(r); r)g-] /4 x 

e r i .rR(Q(r))/12 e -l5(Q(2r),Q(r);r) ^ /jgN 

where R is the scalar curvature of the surface, S the action and D the van Vleck's 
determinant 

d 2 g(Q(2r),Q(r);T) 
' dQ^2r)dQ^r) ' { ' 

£>(Q(2t),Q(t);t) = detail . (21) 

For example for free particles 

TV JV 

« = 2 E ^ = 2 E (<k)« ft , (22) 

t=l i=l 
1 N 

5(Q(2r),Q(r);r) = AT(Q(2r), Q(r); r) = - ^ rf 2 (q,(2r), q,(r))/r , (23) 

i=l 

and for the fluid 

S(Q(2r), Q(r); r) = K(Q(2r), Q(r); r) + ry(Q(r)) . (24) 
We then find the partition function through the integral 

z = J p D (QM;f3)VtidQ, (25) 

and the number density by 

v^(qi)p(qi) = N ^ • (26) 

It is then apparent that by choosing the same isometry on each imaginary time slice we 
reach the same conclusion of Section [3] as for the classical (non-quantum) fluid. 



5. Conclusions 



We showed that in a surface of constant curvature without boundaries the local number 
density p(q) of a non-ideal, (V ^ 0), fluid is a constant in the thermodynamic limit. 
Clearly the ideal gas has constant density on every surface regardless of the curvature 
and of the thermodynamic limit. 

The Coulomb potential for particles living in the surface depends on the metric 
tensor and is in general a function of the geodesic distance between the two charges. 
The Coulomb Quid density is a constant in the thermodynamic limit in the plane [3] 
the sphere [10] (and the pseudosphere [151 ESI E]), surfaces of constant curvatures, but 
also on the Flamm's paraboloid [18], a surface of non-constant curvature. 

We proposed a formulation for the number density which gives to the curvature 
of a surface with a conformal metric (the sphere the pseudosphere and the Flamm's 
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paraboloid are three surfaces of this kind) a physical interpretation as an additional 
external "force" acting on the system of particles moving in the corresponding "flat 
space" . The formulation although suggestive partly masks the intuition of the properties 
of the density because of the fact that the pair potential is inherently related to the 
properties of the curved surface, i. e. the geodesic distance between two points, which 
cannot be translated in terms of the properties of the corresponding fluid moving in the 
"flat space" in a straightforward way. On the other hand the formulation suggests that 
the combination ^fgp is a more fundamental quantity than just p itself. One can show 
both for the non-quantum and the quantum fluid that if d q a is a Killing vector field of the 
surface then if we can neglect surface effects [y/ </(q)p(q)] jQ = and if [\/g(q)p(q)],a = 0, 
V v then g a = 0. These are the main results of our discussion. We can also say that 
g >a = if and only if [\Z^(q)p(q)], a = 0, V v . 

The total potential energy of the fluid moving in the "flat space" is U(Q) = 
V(Q) + J2i j 0) where the functional dependence on Q of the first term depends both 
on the fluid model, through v(d(q^t, cy)), and the kind of surface, through d, whereas the 
functional form of the second term depends only on the kind of surface. It is then to be 
expected that given a fluid model the density can be non-constant on certain surfaces. 

The OCP has uniform density in the cylinder (see Eq. (12a) of Ref. [9]), in the 
pseudosphere (see Sec. 4.3.2 of Ref. [16]), and in the Flamm's paraboloid (see Sec. 4.2.4 
of Ref. [IB] ) • In these cases the explicit expression of the density has been determined 
for the finite system as a function of the properties of the surface at the special value of 
the coupling constant (3q 2 = 2. To the best of our knowledge there aren' t any analytical 
results, in the literature, where the OCP has been found to have a non-constant number 
density in the thermodynamic limit on a given curved surface, probably one has to resort 
to numerical simulations |25j . 

It would be important, in the future, to be able to understand if the surface effects 
on the finite system have some influence in the conclusion that if d q c is a Killing vector 
field of the surface then [-y/ <7(q)p(q)] jQ = in the thermodynamic limit. 
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